Resonance Broadening Induced Nonlinear Saturation of Kinetic Alfven Turbulence in 

the Interplanetary Plasma 
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The saturation of ion cyclotron Alfven turbulence excited by beam particles is investigated 
using resonance broadening theory. The stochastic scattering which decorrelates particles, includes 
both random acceleration by electric fields and a turbulent magnetic mirroring effect. Turbulent 
^ ' mirroring is shown to yield non-Gaussian corrections to the orbits even if the random electric and 

Q magnetic fields are Gaussian. The predicted steady-state turbulence level exhibits a peaked anglular 
■ distribution, with a maximum near O ~ 60°. 
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""q^. It is well known that the interaction of the solar wind plasma with planetary magnetospheres and cometary plasmas 
I ■ results in a high level of wave activity Jjlm[ . Gyrating solar wind ions are a possible agent for exciting high-frequency 
^ instabilities. It has been suggested that these magnetic fluctuations result from ion cyclotron instability, which 
occurs when the solar wind plasma interacts with a beam of cometary ions. A significant effort has been made to 
^ understand the turbulence generated by comet-solar wind interactions using quasilinear theory [4-6] as well as by 
O ^ simulation studies [7[-p^ . Ions are implanted in the solar wind plasma by melting and photoionization, thus forming 
^ I an unstable 'bump-on-tail' like (i.e. shifted Maxwellian) distribution in the solar wind frame ||l^,|lj]. In the far- 
^ r upstream vicinity of a comet, ions form a ring distribution. When the phase velocity of a circularly polarized wave 
coincides with the parallel velocity of particles, i.e. there is an ion cyclotron resonance lo — = Qi between the 
particles and the wave, they resonantly pump the wave. Using quasilinear theory, saturation is predicted to occur via 
plateau formation in the distribution function. This mechanism is rather slow, with a saturation rate of the inverse 
Oh| plateau formation time (i.e. a velocity diffusion time). Saturation via resonance broadening is much faster, with a 
time scale comparable to the instability growth rate (i.e. since the resonance width is proportional to the fluctuation 
^-H [ amplitude). Thus, even if macroscopic modification of the distribution function (a la plateau formation) occurs, it 
^ : should be calculated in the presence of a finite width, fluctuation-broadened resonance. Moreover, the quasilinear 
On ' theory predicts high level of turbulence at ~ 0°, only {Q is the angle between the wave propagation direction and 
the ambient magnetic field), whereas observational data (for example, from the comet P/Grigg-Skjellerup [ p^ , p^ ) 
indicates an the average value of Q between 50° — 60° at large distances from the comet (i? > 400, 000 km. for P/GS) 
and 75° — 90° closer to the comet. The saturation level prediction given here, which is based on resonance broadening 
• theory, yields a peaked 0-profile with the maximum value at Qm — 65° for P/GS parameters. Note that near- normal 
^\ ' angles of propagation can also result from excitation of left hand polarized waves by normal Doppler resonance |l^ ] 
■ in the case when v n c os 6 < va- For Jupiter we predict 0^,^ ~ 60° and observations from the Ulysses spacecraft give 
q; e^P~55°-60° 0. 

' ^ , We develop a nonlinear kinetic model, so it is natural to anzatz that the high-frequency turbulent background is 
composed of a large number of kinetic shear Alfven waves which are excited in plasma by the ion cyclotron instability. 
, The ion-cyclotron resonance becomes w — — fli^ fc|| (va —Vbeam) — ^i. We keep {k±pi)'^ terms, however, since they 
Oh' can be significant at large angles kn/k± <IC 1. An energetic particle immersed in a such turbulent background is subject 

J> . to stochastic scattering which includes both random acceleration by electric fields (E, V|| x B^) as well as a stochastic 
analogue of magnetic mirroring (v^ x B^). Random electric fields result in the diffusion of guiding centers in the 
; ^ ' perpendicular plane. However, the stochastic mirroring effect produces parallel diffusion in velocity space, giving 
rise to superdiffusion and rapid decorrelation (i.e. (a;^) ^ r^) of particle orbit guiding centers along the ambient 
magnetic field. Since turbulent mirroring is an intrinsically nonlinear effect, it results in non-Gaussian behavior orbit 
fluctuation statistics even if the turbulent electric and magnetic fields obey Gaussian statistics. The random Doppler 
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shifts caused by this scattering decorrelate the particles from resonance. To calculate this stochastic decorrelation, we 
employ a modified stochastic orbit theory [|l8l|l9| ] which was developed for both low- frequency (w <C ri^) instabilities 
pO| , pl| and high-frequency (lu ~ fli) instabilities As is shown below, the theory yields results which are in very 
good agreement with the experimental data, in regards amplitude and kj_/k\\ dependence. 

Starting from the Vlasov equation, we separate the distribution function into the averaged (over fast time and 
space scales) part (/(x, v, i)), the test wave response /k(x, v,t), and background responce. We introduce electric and 
magnetic fields of the test wave in terms of scalar and vector potentials as usual, noting _B|| = A± = 0. This 
follows from the fact that kinetic shear Alfven waves are incompressible. Formally applying the procedure developed 
in |2^ ], we now arrive at the lowest order solution of the Vlasov equation 
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where are the factors: /C|| = [l — (1 — fc^w^/tj)(l + k^pD] , /Cj_ = [l ^ + ^IPs)] j ^^^^ Ps ~ Cs/^i, Cs 

is the sound speed, U.^ — eBo/miC is the ion cyclotron frequency. Note that in this formula the resonance function 
R{lu) is defined by the Fourier transform 



/>oo 

R{uj) — / (ir exp(itl}r)(exp [ik • (5x(t 

^0 



r)]) 



(2) 



of the perturbed orbit propagator and asymptotes to the usual resonant denominator of the quasilinear theory when 
fluctuating turbulent forces vanish. Here 5x is the deviation of a particle orbit from the unperturbed one, i.e. 
(5x(t — r) = x'(t — t) — x(i — r) and Q ~ lo — fc||W|| — nrii. Using Eq. (^), we may find (from the dispersion relation) 
the linear growth rate of ion cyclotron instability 
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where no and Ti are the density and temperature of a bulk plasma, and are those of the beam, v/^ is the 
longitudinal beam drift velocity such that ((w — k\\VD — ^i)/k\\VT^^ <C 1, and VTf, is the thermal velocity of the beam 
particles. 

To explore the finite amplitude regime, note that the resonance function R is expressed via the perturbation of the 
orbit function of a test particle, Eq. (y). Expanding (exp(ik • (5x)) in a Taylor series, with the conditions fJ^r 3> 1 and 
T Tkinetic {Tkinetic IS the characteristic time of the quasilinear diffusion), we have 
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(exp(ik • (5x)) = exp{(ik • (5x) - ^((k • A^x)^) + . . .} 
— eyij>{i5uJT — k\D±T — k'uDv..T^ -f- . . .}, 



(4) 



where A(5x — 6x — (<5x), Auj± = k^D±, D± is the perpendicular diffusion coefficient and Dy^^ is the diffusion 
coeflacient of the longitudinal velocity in velocity space. The cubic in r term gives rizc to a F-function in the real- 
space representarion of the resonance function. For the intermediate time-scale ^ t <C Tkinetic as above, one can 
simplify the result by assuming r'^fc^ZJ^n — > T(fc^_Di,|| )^/'^. The resonance frequency shift Suj is proportional to spatial 
derivatives of the guiding center distribution, and vanishes for a Maxwellian |^^. Then, using the approximations 
mentioned above, the resonance function becomes 
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where Au!± = k'^D± and Awy = (fcjjD^n )^/^. 

The system of equations describing the perturbed motion of a test particle cannot be solved exactly because of 
existence of the nonlinear coupling term v x B. We use a perturbation method, assuming small perturbations. Up to 
the second order in perturbation, we have 
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where A^/; is the angle between kj^ and and 

m Jq Jq ml c . 

Here is the unit vector in the plane perpendicular to the ambient magnetic field parallel to vj_, ey is the unit vector 
coUinear to the ambient magnetic field, — ^ iVy and w^q, f ||o are the zeroth order velocities. From here we see 

that even if the £ oc 0k held is Gaussian, J- is not always Gaussian, but rather given by a distribution proportional 
to a Gaussian multiplied by an error function. This means that the third and other moments in the expansion (^) 
are non-zero. The error function is an increasing function of its argument. Thus, corrections appear on the tail of the 
Gaussian distribution, i.e. at sufficiently large times such that r ~ Tkinetic- This effect implies temporal intermittancy 
for strongly mirrored particles. On intermediate time-scales ^ r Tkinetic, this non-Gaussian correction is small 
even when SB / B 1. To capture the simple physical effect of resonance broadening we can drop the non-Gaussian 
terms J- as small perturbations. After straightforward calculation using Eq. (^|J^) and omitting rapidly oscillating 
terms in r, we obtain 

^ [J2_i(s) + 2j2(s) + J2^i(s)] R{u; - /C||«|| - nn,) (8) 



4 



and similarly 
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+ 17^4 [Jn-ii''^) ~ ■^n+M Ri^ - hn - )\ ■ (9) 

4 fc'j,i;^ J J 

The sums in these equations run over the full range of the background turbulence spectrum —oo<k'<oo. Here 
'unprimed' wave vectors refer to the test wave. Eqs. (||j9|) describe the general case of resonance broadening due to the 
effects of scattering of a test particle by random fields. They can be simplified by assuming a narrow-band spectrum 
p2[ . As it will be shown later from estimating the spectrum at saturation, the narrow spectrum approximation is 
valid for this case. Then J^i^ l-S^^I^ /^o^ can be replaced by (5B/B)^. As we consider the ion cyclotron resonance 
with ions we may set n = 1. n ^ I terms give only a small correction. 

The broadening Aa;(v) should be averaged over the velocity distribution function p^ . Saturation occurs when 
the resonance broadening becomes comparable to the linear growth rate 70, i.e. when the total nonlinear growth 
rate vanishes, jnl = 7o ~ = 0. This saturation mechanism does not rely upon quasilinear plateau for- 
mation and thus is applicable to instances where the unstable distribution is maintained externally, such as 
in the case of comet-solar wind interaction. For the ion cyclotron instability, we define the average broaden- 
ing Auj = rj"^(sj) / dvi^M(v)Ji (/i:±-u±/rJi) Aw(v), where = k^vrj^i and Tn{s^) = / d-vFM{^)Jl{k±v±/ni), 
?i = 0, 1, . . .. Upon averaging of Eqs. (||,^ with R{lo) written for the unperturbed case, we can calculate the nonlinear 
saturation level. We use a perturbation analysis to explore those regions where the turbulence level is (sufficiently) 
small, i.e. at angles of propagation close to 0° and 90° (see below). For the case of the perpendicular diffusion (i.e. 
7o — Aw^ =0), we obtain 

5B\^ _ {u:-k\\VD-^^f{l + klpl) ^^^^ 



B 



NL± 1,2 ,,2 



— (1 + fclp )^'^ Fi(sO 



where v\ = Bq/Attp is the Alfven speed, Fi(sj) = ^(s^) / O.SJf (Jq -I- 2J^ + J|) i^M(v)dv. The saturation level is 
roughly proportional to fcj^ ~ sin~^ Q. Thus this case corresponds to the saturation at large angles. By contrast, if 
there is only parallel diffusion (i.e. 70 — Aa;|| =0), we have 



3 




8 [degrees] B [degrees] 

1 Fig. 8 



FIG. 1. The angular distribution of the saturated Alfven turbulence fluctuation level. 
FIG. 2. The angular distribution of the MHD activity level (in arbitrary units) near Jupiter obtained from Ulysses data. 
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where i^2(si) = ^~^{si) J J^FM^dv, ^3(5^) = T-^{si) J 0.25Jf [j§ - j|] FM{v)v±dv. The fluctuation level SB/B 
in this case is proportional to ^ cos~^ &. Thus parallel diffusion suppress turbulence at small angles. In the 

derivation of this equation we made the approximation w f^, but indeed they differ only by a factor close to 
unity. The Awn effect is dominant at small angles of wave propagation, whereas the Au!± effect is more important 
at near-perpendicular angles. For the case where Auj — Auj± + Aa;|| it is necessary to solve the third order algebraic 

equation of the type 70 = ci6B^ + C2&B^^^ . Fig. |l| represent the solution of this equation for typical comet and solar 
wind parameters: n = 10 cm~^,nb = 1.7 cm~^,Te = = 10 eV,Ti, ~ 3 eV,VTi — 3 ■ 10^ cni/s,VA — 0.2vTi,VD = 
2vA, = 10 s-\k = lO-^cm-i. 

Fig. 1 shows the anisotropy of magnetic fluctuation level in propagation angle Q (tanO = k±/k\\). One can see that 
at small angles (0 ~ 0°) with respect to the ambient magnetic field, very small background turbulence levels saturate 
the ion cyclotron instability. The same is true of the near-perpendicular case, 90°. There the turbulence is 
strongly self-suppressed. This is not the case for 8 ~ 60°. Here, the amplitude of the background Alfven waves must 
be large for resonance broadening to self-consistently suppress the ion cyclotron instability. The main energy content 
of the observable MHD activity, therefore, is associated with Alfven waves propagating at large angles relative to the 
averaged magnetic field. The peak is strongly localized at angles near Qm — 60°. The value Qm of the maximum 
turbulent fluctuation is remarkably insensitive to plasma parameters such as bulk plasma and ion beam temperatures, 
density and drift velocity of the beam, etc. For all reasonable parameters, 0^ remains approximately 60 degrees. 
The magnitude of the maximum fluctuation level, on the contrary, depends sensitively on the beam and plasma 
parameters. Fig. ^ is a plot of experimental data of MHD activity near Jupiter, obtained by the spacecraft Ulysses 
p^ . Here the level of fluctuations (in arbitrary units), measured for different angles, is shown. It is apparent that 
the dominant waves are those propagating at angles in the interval 50° < < 70°. There is remarkable agreement 
of the angular distribution obtained from experiments |l^,|l^,Q with our very simple theory. Quasilinear theory 
predicts a different result PJT^], which is (approximately) {•mini/mpnp)'^/'^{v£) /vaY^^ ■ The 0-dependence predicted 
by quasilinear theory is not in agreement with experimental data, although it gives a reasonable prediction of the 
magnitude, i.e. 6B/B ^ 1. 

We emphasize the fact that the peaked angular distribution is the result of both perpendicular and parallel diffusion. 
The perpendicular diffusion (the diffusion of guiding centers) is a result of the magnetic field line flutter and is 
responsible for the resonance broadening at larger angles, Q ^ 90°. The parallel diffusion (the randomized step-size 
of a particle helical trajectory) is due to the random mirroring force. It is responsible for the instability suppression 
at small angles of wave propagation, i.e. for Q near 0°. It is important to note that the parallel diffusion in 
velocity space results in a {x'^} ^ type decorrelation process in real space. However, for intermediate time-scales 
^ T <C Tkinetic, it Can bc modelled by a standard diffusion process. The magnetic mirroring gives rise to 
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non-Gaussian dynamics of a particle immersed into turbulent background, even if this background is Gaussian. This 
describes temporal intermittancy in a system with strongly mirrored particles {SB/ B ^ 1). 
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